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ABSTRACT 

The current standard time delay formula (CSTD) in gravitational lensing and its 
claimed relation to the lens equation through Fermat 's principle (least time principle) 
have been puzzling to the author for some time. We find that the so-called geomet- 
ric path difference term of the CSTD is an error, and it causes a double counting 
of the correct time delay. We examined the deflection angle and the time delay of a 
photon trajectory in the Schwarzschild metric that allows exact perturbative calcu- 
lations in the gravitational parameter GM in two coordinate systems - the standard 
Schwarzschild coordinate system and the isotropic Schwarzschild coordinate system. 
We identify a coordinate dependent term in the time delay which becomes irrelevant 
for the arrival time difference of two images. It deems necessary to sort out unam- 
biguously what is what we measure. We calculate the second order corrections for 
the deflection angle and time delay. The CSTD does generate correct lens equations 
including multiple scattering lens equations under the variations and may be best un- 
derstood as a generating function. It is presently unclear what the significance is. We 
call to reanalyze the existing strong lensing data with time delays. 
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1 INTRODUCTION 

Shapiro (1964) championed the time delay of the radar echo 
off of the inner planets as a fourth test of Einstein general 
relativity. The so-called classic three tests were suggested 
by Einstein and they are Mercury precession, astrometric 
shifts of the bacltground stars due to the gravitational lens- 
ing by the Sun, and gravitational redshifts (Pound and Re- 
bka 1959). In the same year, Rcfsdal (1964b) suggested to 
measure the time delay difference (or arrival time difference) 
between two images to use in conjunction with other mea- 
surables (position and flux of the images) of supernovae to 
determine the Hubble constant and mass of the lens galax- 
ies. The time delay measurement is a central issue in strong 
gravitational lensing where multiple images are detected, es- 
pecially because the time data can help break the degener- 
acy of the lens fitting models. There have been monitoring 
campaigns of multiple image quasars in radio and optical 
(Vuissoz ct al 2008; Coles 2008). With largo high cadence 
survey telescopes in plan (LSST, WFIRST, EUCLID) time 
delay measurements will become a routine business, and it is 
expected to see many multiply imaged supernovae (Kirkby 
et al. 2006). 

The time delay function (2d scalar function) is also 
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known to generate the lens equation (2d vector equation) 
through the stationarity or oxtrcmum condition (commonly 
referred to £is the Fermat's principle) for single scattering 
lenses (Schneider 1985; Read et al 2007) and multiple scat- 
tering lenses withal (Blandford and Narayam 1986). The 
single scattering lens equations are well known from direct 
derivations of the photon path in the Schwarzschild metric 
and its generalizations in the effectively Newtonian gravita- 
tional field. The multiple scattering lens equation can also 
be derived directly by looking at the photon path in 3-spacc 
(Rhie and Bennett 2011) and it confirms the validity of the 
variational method. 

The current standard time delay function is made of 
the geometric path difference term and gravitational poten- 
tial term. The so-called geometric path difference term of 
the CSTD is essentially a quadratic function of the posi- 
tion difference between the source and image, and we have 
been puzzled by it for some time even though it is widely 
used unsuspected and being incorporated into pipeline codes 
for systematic studies of gravitational lensing (Vuissoz et 
al (2008), Coles (2008), Read et al (2007), and references 
therein). We examine the CSTD for a Schwarzschild black 
hole lens which has the virtue of allowing exact calculations 
of the time delays and find the CSTD wrong. It effectively 
doubles the correct time delay. We trace the origin of the 
current standard time delay formula form to Cooke and Kan- 
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towski (1975) (CK75 hereafter) and find that the so-caUed 
geometric path difference term suffers from mistakes. CK75 
uses the so-called isotropic metric (in the linear approxima- 
tion) . Thus we discuss the deflection angle and time delay of 
a photon path in the Schwarzschild metric in two coordinate 
systems - the standard Schwarzschild coordinates and the 
isotropic Schwarzschild coordinates, where the latter can be 
obtained from the former (or vice versa) by changing the 
radial coordinates. In the linear order in rs = 2GM where 
AI is the mass of the Schwarzschild black hole, the radial 
coordinates only differ by a constant: r = u + rs/2 where 
r and u are the standard and isotropic radial coordinates 
respectively. We will see that the two metric forms (both 
asymptotically flat) result in the same time delay difference 
for two images. At the same time, it leaves a question as 
to what is the time the observer measures. We examine the 
deflection angle and flnd that the difference between the co- 
ordinate systems is effectively of the second order. In other 
words, the deflection angles measured in the two coordinate 
systems are the same while the time delays are different. We 
calculate the second order corrections to the deflection angle 
and time delay. We conclude in section 6. 



2 THE CURRENT STANDARD TIME DELAY 
FORMULA 

The current standard time delay formula used in the lensing 
community can be found, for example, in eq.(Al) of Read 
et al (2007) which is reproduced here for convenience. 
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The first term is referred to as the time delay due to the 
geometric path length difference between the deflected and 
undcflected paths, and the second term is referred to as the 
time delay due to the gravitational potential acting on the 
photon along the path. 

For a Schwarzschild black hole of mass M, eq.(l) be- 
comes 
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where Dl is the distance (from the observer) to the lens 
and D is the reduced distance. 9 is the angular position 
variable, and the subscripts s, I/, / denote source, lens, and 
image respectively, zl is the redshift of the lens, (c is the 
speed of light, which will be set to be 1 hereafter: c = 1.) 
We may say that the time function is made of a "quadratic 
term" and a logarithmic term. The stationary condition with 
respect to the variation of the image position, dt/d9 = 0, 
generates the lens equation. Since ^ is a two-dimensional 
angular variable in general, we take it as a vector, and the 
resulting equation is the well-known single lens equation. 
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where 9e is the angular Einstein radius. This stationarity 
condition is commonly referred to as Fermat's principle: an 
image forms such that the time is an extremum (even though 
Fermat's principle for Snell's law invokes the notion of the 




least time.) The single point lens equation (3) has two images 
whose positions are collinear as is well known. 
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9l has been set to zero by translating the coordinate system. 
The arrival time difference {t{9-) — 1{9+): 9- is the position 
of the dimmer image and arrives later) due to the quadratic 
term and the logarithmic term are 

a2 
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Ati2{log) = (1 + zl)8GM In J + (^ + 1)^/2 j _ 

For s < 1, the both become Atimages « (1 + zl)AGMs. 
The total time delay between two images would be ^ (1 + 
zl)SGMs. 



2.1 Refsdal's Time Delay 

Refsdal (1964b) shows that the arrival time difference of 
the two images of a spherical galaxy lens is w 8GMs in 
eq.(6-refsdal). The time delay formula was derived in cq.(30- 
refsdal) of Refsdal (1964a), and wo may write it as follows. 
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where Z3,, = Di + Ds and s is the dimensionlcss source 
position variable. Refsdal didn't present the exact integral, 
but it is easy to do. 
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6+ is nothing but the position of the positive image (in units 
of the Einstein ring radius). The negative image position 
b- = s — b+ and b- = 67^. Thus, 
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It is a surprise; it is exactly the CSTD at zl ~ 0!! Unfor- 
tunately, we have no clue how to understand the definition 
of the time delay in eq.(7): how integrating the angular sep- 
aration of the two images as the source moves from s = 
(where there is no arrival time difference between the two 
images) to s = s is supposed to correspond to the time delay. 
Refsdal doesn't explain it. It must have been self-evident to 
him. In short, we don't know how he got it wrong. 

2.2 Simply Looking at Fig.l 

Figure 1 shows a deflected photon path in (r, (p) scattering 

plane. If we look at fig.l, the path difference between the 
actual photon path and the would-be straight line is A « 
roStj) ~ 2rs. It is a constant irrelevantly of the path. 

We can try to do better. Let's assign viewing angles by 
the observer of the image, source, and lens (M) as a, P, and 
7 respectively, then wc may estimate it as A ~ ri(a — 'y)S(j>. 
Or, we may estimate it a bit more carefully by looking at the 
outgoing and incoming paths separately: Ai w ri(a — 7)^ 
and A2 « r'2(ri/r2)'(Q; — 7)^. Then 
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Figure 1. Depiction of a photon path (arrowed thick curve), from a source 5" to an observer O, passing by a spherically symmetric mass 
M in Schwarzschild coordinates, rg is the distance of the closest approach, Scj) is the scattering angle, and b is the impact distance. In 
the small scattering approximation relevant in most of the astronomical practices, Srji << 1, b a ro and -y/r 
the photon from S to the center of the mass M in the absence of the gravitational interaction. 



Q is the path length of 



A = Ai + A2 
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The path length difference is larger for the brighter image. 
In other words, the arrival time difference of the two images 
will be opposite to what is implied in the "quadratic term." 
The consequence is that the arrival time difference between 
the two images will vanish because the magnitudes of the 
"quadratic term" and logarithmic term are practically equal 
as we have seen in the previous section. The lens equation 
(the relation between the source and image positions) can 
be found easily. 
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Now we can express the geometric path difference A in terms 
of (a-/3). 
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It may be best called an "inverse quadratic term" than 
"quadratic term." 



2.3 Gott and Gunn, Cooke and Kantowski, and 
Schneider's Fermat Principle 

Gott and Gunn (1974) wrote down without explanation the 
quadratic part of Refsdal's time delay in eq. (9) as the cor- 
rect time delay between two images. It is hard to know the 
underlying reasoning. 

Cooke and Kantowski (1975) (CK75 hereafter) clearly 
lays out a reasoning for their calculation and hence is of in- 
terest here. CK75 interpreted the "quadratic term" of Gott 
and Gunn as due to the geometric path difference and pur- 
ported to "complete" the time delay formula leading to the 
CSTD for an arbitrary mass distribution where the weak 
gravity approximation (linear in GAd) holds. CK75 doesn't 



mention that the CSTD was obtained by Refsdal for a spher- 
ically symmetric mass even though they cite Refsdal's 1964a 
and 1964b. It is possible that CK75 didn't know it, and they 
don't mention what might be the underlying reason of Refs- 
dal's definition of the time delay perhaps because of it. CK75 
does offer their own reasoning for the "quadratic term" as 
has been stated already and lays out a rule for the linear cal- 
culation. The "quadratic term" they obtain is of the second 
order in GM, hence it is a violation of their first order ap- 
proximation assumption. Their calculation of the "quadratic 
term" amounts to a confusion where 7 in eq.(lO) is replaced 
by /3, which led to the "quadratic form" instead of the "in- 
verse quadratic form." The difference of factor 2 comes from 
the slight difference in the definition of the path difference, 
(See Fig. 2 of CK75.), which doesn't concern us especially 
because the "quadratic term" will be shown to be irrelevant 
to the time delay. 

Schneider (1985) "improved" CK75's CSTD and intro- 
duced a variational principle to derive the lens equation from 
the CSTD time delay. Schneider called it Fermat's Principle. 
Once the CSTD is found invalid as a measure of the flight 
time of a photon, the notion of Fermat's principle may be 
best abandoned. 



3 TIME DELAY DUE TO A SPHERICALLY 
SYMMETRIC MASS 

3.1 The Standard Schwarzschild Metric 



The deflection angle of a photon trajectory due to a 
spherically symmetric mass M and its flight duration can 
be calculated exactly by solving the equation of motion in 
the following standard Schwarzschild metric. 

ds^ = -B{r)dt''+A{r)dr^+r''dn^; A{r)~'' = B{r) = i_2GM^ 
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The equation of motion can be derived from the variational 
principle, = 5 J ds. 



da" dx^ 
' dp dp 



(14) 



where p is the line parameter and is the Christoffel 
symbol that is a function of the linear derivatives of the 
metric components. 



Kf) = l9''\9^9.xd.g^x - dxg^.) (15) 

where g''^ is the inverse metric components: g^^gxa ~ <5a- 
Because of the spherical symmetry, the photon path lies on 
a plane, say 9 = n/2, and the equations involve three vari- 
ables: t, r, and d>. 
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where J is the angular momentum. From the equations, one 
can get dcjj = g{r)dr or dt = f{r)dr. The scattering (or 
deflection) angle is obtained by integrating d(f) = g{r)dr. 
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The integration J dcj) = J g{r)dr doesn't relent to a nice for- 
mula unlike in Newtonian but can be done easily by expand- 
ing in GM (weak gravity approximation), which is suitable 
for astrophysical purposes. The scattering angle comes out 
to be 5(t> = 4GM/ro as is well known, and the weak grav- 
ity approximation is also commonly called the small angle 
approximation, ro is the distance of the closest approach of 
the photon path to the mass M as is shown in figure 1. In 
order to obtain the the Schwarzschild coordinate time dura- 
tion it takes for the photon to travel from the source to the 
observer, one can integrate dt = f{r)dr in the same small 
angle approximation. 
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where Ba = B{ro). The duration cither from source at r = 
r2 to ro or from ro to the observer at r = ri is given as 
follows. (See, for exmaple, eq.(8.7.4) of Weinberg (1972); the 
speed of light c = 1.) 

t{r,ro) = - rg + 2GMln ^_±V^1EA^ 

+GM{^y\o{rl/r'o). (19) 

where Vs = 2GM is the Schwarzschild radius. The Shapiro 
time delay formula of the radar echo eq.(l-shapiro) of 
Shapiro (1964) can be obtained from eq.(19). 

Note that the first term does not depend on the mass 
M and should be the time the photon takes to travel from 
the source to the center of the mass M when the gravitation 
is turned off. It can be directly calculated from eq.(18) by 
setting A = B = 1 and should also be intuitively evident 
from figure 1. The proper time duration measured by the 



observer is obtained by multiplying the time dilation fac- 
tor B{rY^'^ in eq.(13) to the Schwarzschild coordinate time 
(e.g., eq.(6.3.46) of Wald (1984)). Thus, the total proper 
time delay due to the general relativity in the linear order 
in ra/ro is 
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where r^ and ri are the radial positions of the source and 

observer respectively. Now define the distances from the ob- 
server to the lens and the source along the horizontal line 
D\ and D'i as shown in figure 1. Using 



D]^rl 



where Dj » ro-, 



we get 
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where D is the reduced distance. (Di is used in place of 
Dl in section 2 for handy manipulation of the index.) If 
we ignore the third term for now because it is usually very 
small, the proper time delay between two images, 1 and 2, 
depends only on the first term. 



Ati2 = Ati - At2 = -2rs ln(roi/ro2). 



(23) 



This is exactly the same as the arrival time difference due 
to the logarithmic term of the time function in eq.(2) with 
zl = 0. (Here the expansion of the universe is ignored.) The 
true time delay formula doesn't have a "quadratic term." If 
we recall that A.t{quad) w At{log) when the source is within 
the Einstein radius from the lens (where the dim image is 
not too dim) , the current standard time delay is about twice 
the true time delay. It is not a small difference and it makes 
it urgent to reanalyze the time delay data. Considering that 
the time delay and lens mass modelings have been producing 
reasonable Hubble constants while using a wrong time delay 
formula, it warrants special scrutinies of the fidelity of the 
analyses. 

3.2 The Isotropic Schwarzschild Metric 

The so-called isotropic Schwarzschild metric is of the follow- 
ing form, 

ds^ = -G{u)drf + F{u){du + udQ^), (24) 

and can be obtained from the standard Schwarzschild metric 
by a coordinate transformation: r w. 



u = I (r - r^/2 + (r^ - rr,)^/^) 
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The equations of motion arc 
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The defiection angle and the fiight time of a photon trajec- 
tory are 
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If one examines the deflection angle integral (17) with the 
factor fs in mind, half the deflection angle comes from the 
term A in the numerator and the other half comes from 
the B term in the denominator. F = G~^ in the linear 
order, and one can sec immediately that cq.(28) produces 
the same deflection angle as cq.(17). If we compare the time 
integrals eq.(18) and cq.(29), the latter has the combined 
factor G/F w G^ instead of B in the denominator and the 
extra factor doubles the third term in eq.(19) resulting in 
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In fact, rj{u,uo) can be obtained from t{r,ro) by using 
eq.(25). In the linear order in Vs, u = r — rs/2, and the 
extra factor in the third term in eq.(30) can be seen coming 
from the first term in eq.(19). 
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So we can say that it is algebraically consistent, but it leaves 
a question as to which is the time the observer's clock will 
be measuring. As for the arrival time difference between two 
gravitationally lensed images, the third term is irrelevant in 
the linear order because it is effectively a constant. 
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where the expansion is made in uo/u « 1. 

Incidentally, it is worth noting that ro/r is of the same 
order of Vs/ro for bright images. They are near the Einstein 
ring, and hence r§ ~ R% = 4GMD « 2rsr (see eq.(3)). In 
other words, r^/r — 0{r'^/rQ) for bright images. As for the 
gravitational lensing by the Sun as observed from the Earth, 
the Einstein ring radius is much smaller than the size of the 
Sun {Re ~ -Rq/23.2), and the only visible "image" is highly 
unmagnifled and 1 >> ro/r » Ts/to. 



Note that (r^ 



2x1/2 



in eq.(31) is the distance of the 



straight path the photon would have flown when M = 0. 
So, let's reexamine eq.(17) to see if the angular span of the 
straight path generates an r^-dependent extra term upon 

the coordinate transformation. 
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The first two terms are the angular span of the straight path 
from r to ro, and the second term does generate an extra 
term. In the linear order in rs. 
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It is then curious because the extra term doesn't show up 
in the first order integration of eq.(28). 
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We can reason that it is consistent because the extra term 
is effectively of a second order. 

We will see in a later section that the second order pertur- 
bation in Ta does generate terms of the form r1 / ro which is 
of the order of r^/r for bright images. 

Upon looking at eq.(30), we may reason that the extra 
time is due to that w (u < r) measures a path that is closer 
to the mass where the potential is deeper and so the photon 
speed is smaller. However, it still remains a question which 
formula for an observer to use between cq.(19) and cq.(30) 
to compare with its clock time. (The difference between the 
"observer" 's coordinate time and proper time is the same in 
the linear order in the two coordinate systems as one can 
see in the fourth term of eq.(20).) Given a mass M, which is 
rightfully considered a coordinate- invariant, wc can imagine 
an abstract space and two abstract points S and O on it and 
a null geodesic connecting the two points. A coordinate sys- 
tem is a means to describe the abstract situation. What we 
are at a loss seems to be whether we are identifying the same 
abstract points S and O in the two well-known coordinate 
systems in which we examined the time delays and deflection 
angles. It should be worth sorting it out unambiguously. 



4 GEOMETRIC PATH DIFFERENCE TERM 
AND GRAVITATIONAL POTENTIAL TERM 

Let's write the isotropic Schwarzschild metric in eq.(24) in 
the linear approximation as CK75 did for more general grav- 
itational fields. 

ds" = -(1 - -)drf + {1 + -)(du + udif). (36) 

M W 

(CK75 uses (two times) the Newtonian potential — 2$(u) 
instead of rs/u.) For a photon trajectory, ds^ = 0, and the 
time integral can be written formally as follows. 



dri = dss + 



-dsi, 



(37) 



where dss denotes the length element of the photon trajec- 
tory in 3-space (to distinguish from the four dimensional line 
element ds). CK75 interpreted, 

"The first term is the time due to the length of the path 
traveled and must be computed to first order in G 
[our highlight]. The second is due to the potential well 
through which the photon traveled." 

It is intuitive and reasonable. As far as we know, this is the 
origin of the current standard notion of the time delay being 
constituted of two pieces: path difference and gravitational 
potential. 

The question is what dss is, and we need to solve for 
the trajectory to get the accurate coordinate time delay. 
It is instructive to consider the straight path. If dss is the 
line length of the straight path, then J dss = (r^ — ^o)^^^ 
and hence ds3 = (r^ — ro)~^^^rdr; the integral eq.(37) re- 
produces the first two terms of eq.(30). The third term of 
eq.(30), which we discussed to be coordinate dependent and 
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irrelevant to the time delay between two images, is not re- 
covered. For the Schwarzschild metric, we we have solved 
the equations of motion and know that the correct line el- 
ement is dp. It reproduces the time integral eq.(29 exactly 
in the linear order. The integral eq.(37) can be rewritten as 
follows in the linear order. 

where £ is the line parameter of the straight path. We can 
see that J dp generates the geometric path difference effect 
(« 2rs) we have seen in section 2.2 and CK75 might have 
liked to have. 

As we have repeatedly stated, the path difference term 
2rs cancels out for the arrival time difference of two images 
for a Schwarzchild lens. Thus for a general potential where 
the form of the correct line element dp may not be easily 
obtained, it may be reasonable to use the straight line ele- 
ment. 

In section 2.2, we "tried to do better" and obtained the 
"inverse quadratic term" which we know to cancel out the 
true time delay. It is not any better than CK75's doubling. 

As we have been repeating it, the magnitude of the "inverse 
quadratic term" is approximately 2rs ■ The fault seems to be 
lie in that we introduced the second order in angles while 
the premise is the linear order small angle approximation. 
CK75 erred the same. The geometric relation rs/ro ~ ro/r, 
where r is the position of the source or observer, seems to be 
the source of ready confusion. The lesson may be that one 
needs to be as systematic as possible. There is no obvious 
way to obtain the correct line element dp without solving 
the equations of motion, and the best bet is the straight 
line element. It is unlikely that CK75 would have been able 
to get the correct geometric path difference. However, the 
notion of geometric time difference is fine. 



5 THE SECOND ORDER CORRECTIONS 

The second order terms of the flight time and angular span 
are presented here. If we define 
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the part of the time integral eq.(18) that is of the second 
order in is 
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If we set r = oo, At2{r = oo,r()) — 2.445 r'i /ro. It is good to 
find the coefficient of the second order to be of 0(1). The 
second order part of the deflection angle eq.(17) is 



6 CONCLUSION 

We examined the current standard time delay formula for 
a Schwarzschild black hole lens in the weak field (or small 
angle) approximation. We find that the current standard 
time delay formula is wrong. It effectively doubles the true 
time delay. The 'quadratic term" is the result of a wrong 
calculation of the geometric path difference. 

The absence of the "quadratic term" thwarts the claim 
of the relation between the time delay and the lens equation 
through the Fermat's principle. However, the CSTD formula 
can be considered a generating function of the lens equations 
because it works, even though it is unclear presently what 
significance it has. 

We call to reanalyze the strong gravitational lenses with 
time delay measurements. We expect the correct time delay 
formula to help for better lens modelings. On the other hand, 
the current practices of reasonably good fittings to the Hub- 
ble constant while using a wrong time delay formula makes 
one worried lest the data analyses tend to converge on desir- 
able values. We call for exceptional scrutinies of the fidelity 
of the analysis methods. 

We identified a term in time delay that is dependent 
on the coordinate systems, even though it is benign for the 
arrival time difference of two images. The coordinate depen- 
dence of the "measurables" seems to require some careful 
thoughts as to what is what we measure. 
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